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Abstract 

The deuteron magnetic and quadrupole moments are unambiguosly determined 
within the front-form Hamiltonian dynamics, by using a new current operator which 
fulfills Poincare, parity and time reversal covariance, together with hermiticity and the 
continuity equation. For both quantities the usual disagreement between theoretical 
^ ■ and experimental results is largely removed. 
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The deuteron is a good system for a test of relativistic approaches devoted to the investiga- 
tion of hadron electromagnetic (em) properties (see, e.g., |[], |2[ ||, ||, || and Refs. quoted therein) 
and in particular of the accuracy of the one-body impulse approximation (IA) for the current 
operator. It is usually believed that effects beyond IA, e.g. meson-exchange currents, iVA-pair 
creation terms (Z graphs), and isobar configurations in the deuteron wave function are important 
for the explanation of existing data. However, these effects are essentially model dependent || and, 
furthermore, obviously depend on the reference frame (see, e.g., Refs. |?], ||). 

A widely adopted framework for relativistic investigations of deuteron em properties is the 
front-form Hamiltonian dynamics (FFHD) 0, |10| , where only the two-nucleon state is considered 
and the wavefunction of the system factorizes, for any front-form boost, in an eigenfunction of the 
total momentum times an intrinsic wavefunction. Thus, deuteron em form factors, determined from 
three independent matrix elements of the current, are given in terms of elastic em nucleon form 
factors and deuteron internal wave function. In the FFHD, the one-body approximation was usually 
applied only to the relevant matrix elements of the plus component of the current in the reference 
frame where q + = (q is the momentum transfer), while the other ones were properly defined 
in order to fulfill all the required properties (see, e.g., |l], 0, 0). In this letter, we will consider 
the Breit reference frame where the three- momentum transfer is along the spin quantization axis, 
which allows one to exploit the symmetry of the problem and to calculate all the non-vanishing 
matrix elements of the current by the same rules. 

Following Ref.|8|, let us consider the current which in the Breit frame where P± = q± = 
has the form 

f(Ke z ) = ±{J»{Ke z ) + L£[r x (-ir)] ex P (nrS x )J v (Ke z )*ex P (-mS x )} 

J + (Ke z ) = J-{Ke z ) = (P ± = 0,P' + \UJ+ ree (0)U\P ± = 0,P+) 

J±(Ke z ) = (Pi = o,p'+|nj /ree± (o)n|p ± = o,p+>. (i) 

In Eq. (jl|), n is the projector onto the subspace of deuteron bound states |xi) of mass 
rrid and spin 1, J^ ree (0) = J^(0) + J^(0) is the one-body current, |Pj_,P + ) is an eigenstate of 
the total deuteron momentum, P ± = (P x ,P y ), P + = (P + P z )/V^ = -75 [("^ + K 2 ) 1 / 2 - K], 
P'+ = -L[( m J + K 2 ) 1 ' 2 + K], K = Q/2, Q 2 = -q 2 and q = P' - P; L[r x (-n)] is the element 
of the Lorentz group corresponding to a rotation of —it around the x axis, S x is the x component 
of the front-form spin operator, and * means Hermitian conjugation in internal space. From Eq. 
(|l|), one can obtain the expression of the current in any other reference frame by applying the 
proper transformations (see, e.g., Q). This current operator fulfills extended Poincare covariance, 
hermiticity and the charge normalization, as well as current conservation ||. The second term 
in the first line of Eq. (|l|), which ensures hermiticity, introduces two-body terms in the current, 
because of the presence of S x (see below). 

A relevant result of our approach is that the extraction of elastic em form factors is no more 
plagued by the ambiguities, related to the so called "angular condition" , which are present when the 
free current is used in the q + = frame (see, e.g., [II], |7|, fL^l). In this case, one has four independent 
matrix elements of the current, while the em form factors are three [11]. Differently, in our model 
(Eq. (|l])), it turns out || that only three matrix elements jg, s = (mdlS'Jj^^e^lmdlSz) are 
independent (e.g., j$ , j^, and jf ). Therefore, there is no longer any freedom in the construction 
of the three em form factors. These matrix elements (as well as any other one) can be easily 
obtained by Eq. (]!]) in terms of the matrix elements Ja, s = {mjlS'^J^ (K e z )\md\S ' z ) . Indeed, 



by using the properties of the Wigner D-functions, one can show that the two terms in the first 
line of Eq. (||) are equal for the plus component (Jq m ]\ = i7o + 0(l i))> w h ne for the x component 
they yield Jft/2 and -Jfa/2, respectively (jf fi = [Jfr- J^}/2) §. 

As a test of our current, we evaluate the deuteron form factors at Q 2 = 0, namely the 
magnetic moment, Hd, and the quadrupole moment, Qd, which are not affected by the uncertainties 
in the knowledge of the neutron em form factors at finite momentum transfers. The deuteron 
moments are a longstanding problem in nuclear physics, since it was not possible to reconcile in 
a coherent approach theoretical and experimental values for both quantities at the same time, by 
changing the tensor content of the nucleon-nucleon (JV — N) interaction, or considering two-body 
current contributions, both in non-relativistic and in relativistic frameworks 0, 13, O, p|. Our 



preliminary results for the deuteron form factors at Q 2 ^ can be found in [15 1 . 

By using the properties || of the matrix elements of jg, s the deuteron form factors can 

be written in terms of the matrix elements Jgi a [15]. Then, the magnetic moment, in nuclear 

magnetons, is given by 

(m p \/2) 1 [Jiq - Joji] /rA 

Vd = — lim — — , (2) 



while the quadrupole moment is 



V2 .. 1 

hm —, 

m d Q^o Q 1 



g d = £-hm- f [J + - l 7+]. (3) 



If one adopts the free-current in the q + = frame, the angular condition is satisfied at 
the first order in Q, but it is violated at the second order, for Q 2 — ► nj. Therefore the angular 
condition is not a problem for the calculation of jidi while the quadrupole moment is not uniquely 
determined. Differently, following our model, both of them are well determined. 

The matrix elements J% q can be easily calculated, by using the action of the free current 

on a two-body state \P±,P + )\xs,S z ) 0: 

(A^AA\J}ree(0)\P± = 0, P + )\XS,S Z ) = E«(piVl)IVC/? (W ~ Pi) 3 ) " 

0"l" 

mti-Pif ft 1 ±4w + /£((Pi -Pi) 2 )7XPl^l")(fe^l" ; 4lX5A)^ (4) 

where w(p,a) is the front-form Dirac spinor [ff3], while /* s and /^ are the isoscalar electric and 
magnetic Sachs form factors of the nucleon. The relations between the internal (k±, k z ) and 
individual nucleon variables in our reference frame are given by 

Pi±=0 1± = k±, pt=£P + , k g = u(k)(2£-l), £ = 1 + (Z-1)P + /P' + (5) 



where co(k) = \Jrri 2 + k 2 , with m = (m p + m n )/2 the nucleon mass, and k = \k\. Nucleon form 
factors cannot be factorized in the current matrix elements, since from Eq. (||) one has (p'i —pi) 2 = 
-4Q 2 (m 2 + k 2 ± )/(4m 2 tt'). 

In FFHD, the internal deuteron wave function with polarization vector es is given by (cf. 

§) 

1 "^kk- 

<*|X1A> = (27rf 2 Hk)/2} 1 / 2 v(kr 1 v(-k)- 1 [Mk)S l j - ^fo " -^r L )^{k)]a i a y (e Sz ) J (6) 



where a sum over the repeated indices i,j = 1,2,3 is assumed, v(k) is the generalized Melosh 
matrix ||] and <r, are the Pauli matrices. The wave functions tpo(k) and (f2(k) coincide with 
the non-relativistic S and D waves in momentum representation [|lj] and are normalized so that 
f[cp (k) 2 + ¥2(k) 2 ]d 3 k = l. 

Our FFHD results corresponding to different N — N interactions are compared in Table [l| 
with the non-relativistic ones (for overcoming numerical instabilities a careful analytical reduction 
of Eqs. (|2],||) is needed). The standard non-relativistic results obtained with a one-body current 
crucially depend on the asymptotic normalization ratio rj of D and S wave functions and on the 
D— state percentage, Pd, but one cannot obtain at the same time the experimental values for both 
Hd and Qd- In our Poincare covariant calculation the relativistic corrections (RC) bring both fid and 
Qd closer to the experimental values, except for the charge-dependent Bonn interaction. In Ref. 
0|, RC have been calculated within FFHD by using the free current in the q + = frame and they 
resulted to be very small for Qd, while for fid were able to explain only part of the disagreement 
with the experimental value. It should be stressed that our current operator and the one used in 
Ref. are different, since both of them are obtained from the free one, but in different reference 
frames, related by an interaction dependent rotation. 

In Fig. 1, fj>d and Qd are reported against the asymptotic normalization ratio, r\. As already 
observed for the non-relativistic calculations of Qd [EoL |26]j , a remarkable linear behaviour appears 
for both quantities, except for the Bonn interaction. The values of fid and Qd, suggested by this 
linear behaviour in correspondence of rj exp = 0.0256, differ from the experimental ones only by 
0.5% and 2%, respectively, i.e. much less than for the non-relativistic results. The RC to [id are 
rather large and the total result is greater than fi d xp - This shows that, within our framework, 
even the sign of explicit contributions of two-body currents is different from the one needed in the 
non-relativistic case. 

In summary, our results for fid and Qd, unambiguosly calculated by a Poincare covariant 
current built up from the one-body current in the Breit reference frame where P± = q± = 0, 
show that the total contribution of explicit two-body currents (from meson-exchange, Z-graphs, 
etc.) and isobar configurations is relatively small at Q 2 = 0. It should be stressed that explicit 
two-body current contributions, considered in addition to the ones already contained in Eq.(^), 
must fulfill separately the constraints imposed by the extended Poincare covariance, hermiticity 
and current conservation jj|. An evaluation in our Breit frame of explicit two-body contributions 
will be performed elsewhere. 

The authors wish to thank A. Kievsky for kindly providing the deuteron wavefunctions for 
RSC, Avl4, and Avl8 interactions and R. Machleidt for the CD-Bonn wavefunction. 
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Table 1: Magnetic moment (in nuclear magnetons) and quadrupole moment for the deuteron, 
corresponding to different N — N interactions; fi^ R and Q^ R are the nonrelativistic results, 
Hd, (LPS) and Qa (LPS) our present results; Pp is the .D-state percentage, and r\ = A D /A S 
the asymptotic normalization ratio. 



Interaction 


Pd 


V 


..NH 


Hd (LPS) 


02" Um') 


Q d (LPS) (fm 2 ) 


Exp 




0.0256(4) |8| 




0.857406(1) U| 




0.2859(3) m 


RSC f~ 


m 


6.47 


0.0262 


0.8429 


0.8611 


0.2796 


0.2852 


Avl4 


m 


6.08 


0.0265 


0.8451 


0.8608 


0.2860 


0.2907 


Paris 


23) 


5.77 


0.0261 


0.8469 


0.8632 


0.2793 


0.2841 


Avl8 


m 


5.76 


0.0250 


0.8470 


0.8635 


0.2696 


0.2744 


Nijm93 pi 


5.75 


0.0252 


0.8470 


0.8629 


0.2706 


0.2750 


RSC93 || 


5.70 


0.0251 


0.8473 


0.8637 


0.2703 


0.2750 


Nijml p|f 


5.66 


0.0253 


0.8475 


0.8622 


0.2719 


0.2758 


CD-Bonn [§|] 


4.83 


0.0255 


0.8523 


0.8670 


0.2696 


0.2729 



Figure Caption 
Fig. 1. (a) Deuteron magnetic moment, /x^, as a function of the asymptotic normalization 
ratio r/, for different N — N interactions. The full dot represents the experimental values for /i^ and 
77; empty triangles and diamonds correspond to the non-relativistic and relativistic results of Table 
|], respectively, while the solid and dashed lines are linear fits for these results. Full triangle and 
diamond are the results of the CD-Bonn interaction, (b) The same as in (a), but for the deuteron 
quadrupole moment, Qd- 
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